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Abstract. In this paper we derive the Markowitz-optimal, deter- 
ministic-execution trajectory for a trader who wishes to buy or sell 
a large position of a share which evolves as a geometric Brown- 
ian motion in contrast to the arithmetic model which prevails in 
the existing literature. Our calculations include a general tem- 
porary impact, rather than a specific function. Additionally, we 
point out — under our setting — what are the necessary ingredients 
to tackle the problem with adaptive execution trajectories. We pro- 
vide a couple of examples which illustrate the results. We would 
like to stress the fact that in this paper we use understandable 
user-friendly techniques. 



1. Introduction 

The problem of optimal execution is a very general problem in which 
a trader who wishes to buy or sell a large position K of a given asset 
S — for instance wheat, shares, derivatives, etc. — is confronted with the 
dilemma of executing slowly or as quick as possible. In the first case 
he/she would be exposed to volatility, and in the second to the laws 
of offer and demand. Thus the trader most hedge between the market 
impact (due to his trade) and the volatility (due to the market). 

The key ingredients to study this optimization problem are: (1) The 
modeling of the asset — which is typically modeled as a geometric Brow- 
nian process. (2) The modeling of the so-called market impact which 
heuristically suggests the existence of an instantaneous impact — so- 
called temporary — and a cumulative component referred to as perma- 
nent. And finally (3) one should establish a criteria of optimality. 

The main aim of this paper is to study and characterize the so- 
called Markowitz-optimal open-loop execution trajectory, in terms of 



nonlinear second order, ordinary differential equations (Theorems 3.2 
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and 3.6). The above is done for a trader who wishes to buy or sell 
a large position K of shares S which evolve as a geometric Brownian 
motion (although in the existing literature it is considered an arith- 
metic Brownian motion for this problem). In this paper we only deal 
with deterministic strategies, also called open loop controls; however, 
we point out — in Section [5] — the key ingredients to address the prob- 
lem with adaptive strategies, also termed Markovian controls (work in 
progress) . 

The main motivation of this work is on one hand economic, but on 
the other the effect of market impact in the valuation of contingent 
claims, and its connection with optimal execution of derivatives. Intu- 
itively, for this kind of problem, one would expect to consider adaptive 
strategies to tackle the questions, although it seems natural to first un- 
derstand the deterministic case. An important element in our analysis, 
will be the use of a linear stochastic differential equation, first used — 
to our knowledge — by Brennan and Schwartz (1980) in their study of 
interest rates. 

The problem of minimizing expected overall liquidity costs has been 
analyzed using different market models by Bertsimas and Lo (1998), 
Obizhaeva and Wang, and Alfonsi et al. (2007a,2007b), just to men- 
tion a few. However, these approaches miss the volatility risk associ- 
ated with time delay. Instead, Almgren and Chriss (1999,2000), sug- 
gested studying and solving a mean-variance optimization for sales rev- 
enues in the class of deterministic strategies. Further, on Almgren and 
Lorenz (2007) allowed for intertemporal updating and proved that this 
can strictly improve the mean-variance performance. Nevertheless, in 
Schied and Schoneborn (2007), the authors study the original problem 
of expected utility maximization with CARA utility functions. Their 
main result states that for CARA investors there is surprisingly no 
added utility from allowing for intertemporal updating of strategies. 
Finally, we mention that the Hamilton- Jacobi-Bellman approach has 
also recently been studied in Forsyth (2009). 

Our paper is organized as follows: In Section 2 we introduce our 
model, assumptions and auxiliary results. Namely through a couple 
of subsequent Propositions we characterize and compute — by use of a 
Brennan-Schwartz type process ([s]) — the moments of certain random 
variable that is relevant in our optimization problem. Section 3 is 
devoted to deriving and proving the characterization of our optimal 
trading strategies and the optimal value function as well. After that, 
in Section 4, we present a couple of examples in order to exemplify 
the procedure derived in Section 3. We first compare Almgren and 
Chriss trajectory with ours, and in Example 2 we use a temporary 
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market impact h to the power 3/5 as suggested in the empirical study 
of Almgren et. al.^s (2005). We conclude the paper, in Section 5, by 
pointing out the key ingredients in the study of adaptive execution 
strategies. 

2. Auxiliary Results 

In this section, we describe the dynamics of the asset S, the so- 
called market impact, and we introduce the Brennan-Schwartz process. 
This process will allow us not only to compute the moments of the 
optimization argument, but also to represent it in terms of an SDE. 
For the remainder of this section, let c{t) be a fixed and differentiable 
function for < t < T < cxd. 

2.1. The model. Let the price of the share S" of a given company 
evolve as a geometric process, where the random component B is stan- 
dard Brownian motion, i.e. 



dSt = Si 



g{c{t)) + ^^^^^]dt + adB 



So 



and where g and h represent respectively the permanent — which accu- 
mulates over time — and instantaneous temporary impact. Thus, the 
future effective price per share due to our trade can be modelled as 

(1) S't = sexp< / g{c{v))dv + h{c{u)) — -a'^u + aBu 



where a > is an estimable parameter. 



Remark 2.1. (a) Note that that process In(S'), where S is as in 
coincides precisely with the "standard" notion of both permanent and 
temporary impact. That is, if we model the price process as arithmetic 
Brownian motion: 

dh 

dSt = adBt + g{c{t))dt + —{c{t)), So = s 

then 

St-s= [ g{c{u))du + h{c{t)) + aBt. 
Jo 

Hence, the first term in the right-hand side of the equality is accumu- 
lating over time, on the other hand, the second term is not. 
(b) Next, observe that the process c can be thought of as a control, which 
in turn may be: 
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(1) an admissible process c which is adapted to the natural filtration 

of the associated process ^ is called a feedback control, 

(2) an admissible process c which can be written in the form Ct = 
u{t, St) for some measurable map u is called Markovian control, 
notice that any Markovian control is a feedback control, 

(3) a deterministic processes of the family of admissible controls are 
called open loop controls. 

(c) In this paper we will only deal with open loop controls, yet to study 
feedback controls it will become quite useful to introduce the so-called 
Brennan- Schwartz process, introduced in the next subsection. The rea- 
son being, that to derive the Hamilton-Jacobi-Bellman equation — see 
Section 5 — it is more convenient to have a diffusion instead of an av- 
erage of a diffusion. 

2.2. Averaged geometric and Brennan-Schwartz processes. In 

order to study our problem we will introduce the following averaged 
geometric Brownian process: 

(2) it := / c{u)Sudu, t>0. 

Jo 

By ([T]) we can express ^t as 
io 

Thus, if c{u) represents the number of shares bought or sold at time u 
at price Su, then ^t represents the total amount spent or earned by the 
trader up to time t. 

To compute the moments of we will use the following linear non- 
homogeneous stochastic differential equation 



(3) dXt = 

Xn = 



c{t) - ( 9{c{t)) + -a']Xt 



dt — aXtdBt 



which has been used for instance by Brennan and Schwartz (1980) in 
the modeling of interest rates, by Kawaguchi and Morimoto (2007) 
in environmental economics, and which may also be used to study 
the density of averaged geometric Brownian motion [see for instance, 
Linetsky (2004)]. 

In general, its usefulness is due to the fact that one may construct a 
Brennan-Schwartz process X which satisfies 
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where = stands for equality in distribution. In this paper, it is used, 
together with Ito's lemma to show that = S ■ X, which alternatively 
will allow us to compute the second moment of ^ in terms of an iterated 
integral, indeed: 

Proposition 2.2. Let processes S, ^, and X be as in and 

respectively, then 

^t = SfXt, t > 

and 

o?^2 ^ 2^tStc{t)dt. 

Proof. By Ito's lemma 
diSfXt) = StdXt + XtdSt + dXtdSt 

= Stc{t)dt - StX, (^g{c{t)) + ^^^^^ dt + StX.a^dt 

-aS,X,dB, + XA (^g{c{t)) + dt 

+aXtStdBt - a^StXtdt 
= c(t)Stdt 
= 6- 

Furthermore, for the second moment of ^ it follows that 

d^^ = 2XtSfdXt + 2X^StdSt + Sf {dXtf + X2 {dStf 
+AXtSt{dXt ■ dSt) 

= 2XtS^sc{t)dt - 2XtS^ (g{c{t)) + Xtdt 



dt 

+2XtS^a^Xtdt - 2XtS'^aXtdBt 

+2XfSt (^g{c{t)) + ^'^^^ Stdt + 2X^StcxStdBt 

+Sfa'^Xfdt + X'^a^S'fdt - AXtSta^XtStdt 
= 2XtSfc{t)dt 
= 2^tStc{t)dt, 

as claimed. □ 

Remark 2.3. The previous proposition may be derived directly from the 
integration by parts formula. Yet, this characterization will be useful 
in the study, for instance, of the optimal trading schedule of derivatives 
or in the determination of Markovian controls. 
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2.3. Moments of ^. Now, by Proposition 2^, it is straightforward to 
compute the first two moments of C,t which will be used to solve our 
optimal execution problem. 

Corollary 2.4. Let ^ be as in pi). Then 



(4) m] 
(5) 

Proof. By Q: 

mt] 



c[ujsexp 
t 



g{c{v))dv + h{c{u)) J> du 

c{u)se-^o g{c{n))dn+h{c{u)) 
X ( [ c{v)se-^o9iciw))dw+h{c{v))+a^v^^\ 



c{u)sK[Su]du 



c[u)sexp 



g{c{v))dv + h{c{u)) \ du 



From Proposition 2.2 



2E 
2E 



c{u)sSu^udu 



c{u)s'^Su ( / c{v)sSydv I du 



CMS 







c{v)s'E[SuSy]dv I du. 



Therefore, since 



]g |^gO--Bu-|-cr_B„ j _ ^^a{Bn-By)+2ijBy~^ 



62 



it follows that 



(.(^y^jg^So 9i.c{n))dn+h{c{u)) 

X ( [ c{v)se^o^«^'>^<i^+h«^))+-'^dv ) du. 



□ 
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3. Markowitz Optimal open-loop Trading trajectory 

In this section we derive a Markowitz-optimal open-loop trading 
strategy, Theorem |3.2| and |3.6[ employing the auxiliary results derived 
in the previous section. 

3.1. Execution shortfall. If the size of the trade K is "relatively" 
small we would expect the market impact to be negligible, that is, 
the trader should execute K immediately. Thus, it seems natural to 
compare the actual total gains (losses) E,t with the impact-free quantity 
Ks by introducing the so-called execution shortfall Y defined as 

(6) Y:=iT- Ks. 

If we use Markowitz optimization criterion, then our problem is 
equivalent to finding the trading trajectory {c(i:)|0 < t < T} which 
minimizes simultaneously the expected shortfall given a fixed risk- 
aversion level A characterized by the volatility of Y: 

J[c(-)] := E[F] + AV[F] 

(7) = \neT]+niT]-\mT\f -Ks. 

In fact, if A > then ([T]) has a unique solution, which may be repre- 
sented in the following integral form: 

Proposition 3.1. Suppose that the permanent impact g is linear, i.e. 

g{x) = ax, 

for some a > as suggested by Almgren et. al. (2005) empirical study. 
Let 



(8) fix) ■= / ciu)du, fix) := cix) 

Jo 

and 



T 



XKs 



Jo 

Then J[c(-)] m can be expressed as: 

[2A (7i(n) - -f{u)) + 1] /'(M)se"^(")+'^(^'("» - ^ |> du. 
Proof. Setting 

fix) := I ciu)du, 
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we have f'{x) := c{x). Hence, using the integration by parts formula 







t / 



imphes 



t \ 2 



= mt]f- 

Thus, by Q and Q, 



Jo 

It follows that 

Jo io 



(2A7i(m) + 1 - 2A7(m)) /'(M)se"^(")+^(^'(")) - ^ Um. 

as claimed. □ 

Observe that this last expression has the following functional form 
in terms of /: 

(9) = f cMu,f,f'),^{u,f,r),f{u),nu))du. 

Jo 

Letting 

F{f{u), f{u)) := sf'iu) exp {af{u) + h{f{u))} , 
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we may re-express ^ as 

J(/)= r (2\ r F{f{v)J\v)){e'^'^ ~l)dv + l\F{f{u)J\u))du. 



In particular 

Theorem 3.2. Suppose that A = 0. Then, the open-loop trading sched- 
ule c* is determined by a function fi which solves the following system 

dF d OF 

With /i(0) = and fi{T) = K. 

Proof. If A = 0, that is, if we only wish to minimize expected execution 
shortfall, then: 

J(/)= r F{f{u)J'iu))du, 
Jo 



and thus equation (10) follows from the Euler-Lagrange equation [see 
for instance Gelfand and Fomin (2000)]. □ 

Example 3.3. Let T = a = 1 and both temporary and permanent 
impact be linear, i.e. g{x) = x, h{x) = x, hence: 

F{f{u), f\u)) = f\u) exp {/(«) + /'(«)} . 

Thus, if one wishes to minimize the expected execution shortfall one 
should execute according to: 

nu) - r{u) - (1 + f\u)){f{u) + r{u)) = o, 

given that /(O) = 0, and /(I) = K. 

Remark 3.4. Note that as A increases the client is willing to be more 
exposed to risk. This idea is equivalent to saying that he/she is willing 
to increase the speed of execution. Hence we have that for A > 0, the 
optimal trajectory will dominate fi: 

f{s)>fi{s), 0<s<T 

In other words we may decompose f{s) = + /2(s). This last 

observation together with our constraint /(O) = 0, /(T) = K, lead to 
the following two facts: 

(11) /2(t) > 0, 0<t<T 

(12) /2(0) = /2(T) = 0. 



Remark 3.5. The previous remark and equation suggest a 2 step 
procedure to find the optimal trajectory f . Namely, first find fi and 
given that information solve for / = /i + /2 
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Theorem 3.6. The optimal dijferentiable trajectory f which solves 
is given by fi + f2, where fi is given in Theorem 3.3 and f2 satisfies 
for < V < T : 

dF d fdF 
Iff'd^ [df' 



/2(«)(^2A F\f{v)J\v))dv + l 
+2A / h{v) 



df dv V df 



dvFifiu)J'iu)) = 0, 



where /2(0) = /2(T) = 0. 



Proof. The idea is to foUow the derivation of the Euler-Lagrange equa- 
tion [see for instance, Gelfand and Fomin (2000)], but in this case, the 
unknown function /2 will play the role of the perturbation. Thus, it is 
essential the fact that /2(0) = f2(T) = 0. 
Let 

f{v) = f,{v) + f2{v) 

g,{v) = fi{v) + ef2{v) 
where /2(0) = = /2(T) and 



J(e) 



T 



2A / F{g,{v),g',iv))ie'^^-l)dv + l]Fig,iu),g',iu))du 



then 



dJ^ 
de 



2A 



" dF 
de 



{9e{v),g:{vme^''''~l)dv 



xF{g,{u),g'^{u))du 
+ [ (2X r F{g,{v),g[{vW''^ -l)dv + l 



dF 

x-^i9e{u),g'^{u))du 



2A 

xF{g,{u),g'^{u))du 
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if we set 



7eW = 2A / F{g,{v),g',{v)){e'''''' - l)dv + l 
Jo 

and by the integration by parts formula we have that 



dJ_ 
de 



£ (^2\ iy,{v){e'^'^ - 1) 
-2A / f,(v)—a'e'''^dv 

Jo dg'e 



dg, dv\dg'^ 



dv 



+2XMu) — {e'^'^ - 1) ]F{g,{u),g',{u))du 



+ 



OF d (OF 
% du \dg'^ 



I %{u)f2{u 

Jo 



du 



Let us compute 
d 

du 



{^,{u))^2XF{g,{ulg[{u)W'^-l) 



which yields 

J'{e) = f{2\jy,{v){e^'^-l) 



dF _ d (OF 
% dv\dg'^ 



dv 



f,{v)—a'e'^'^dv j F{g,iu),g'^{u))du 



I 



+ / 7.(«)/2(«) 



'dF _ d /OF 
% du \dg'^ 

But observe that when J'(l) we have 



du. 



J'(l) = ^ i2xjj,{vW'^-l) 



OF _ d (dFY 



dv 



-2^ / f2{v)^a'e^'^dv ]F{f{u),f{u))du 



+ 

= 0. 



/ ll{u)f2{u) 

Jo 



'dF 








du 





du 
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Now, given that /2(0) = f2{T) = and from the fundamental lemma 
of Calculus of variations we have that: 



f2{uhl{u 

+2A 



'dF 








du 





f2{v) 



d fdF^ 



dvF{f{u)J'{u)) = Q 



'OF^ _ 

1)1 ~ dv\ df 

with the constraint that /2(m) > for m G (0, T) and f2{T) = /2(0) = 0. 



or 



Mu)i2\ / F\f{v),f\v))dv + 1 



dF d fdF 



+2 A / h{v) 



'dF^ 


d 


rdF'\ 


[df 


' d^ 


[df'J 



df du \df\ 
dvFifiu)J'iu)) = 



as claimed. 



□ 



Example 3.7. (cont.) With linear and temporary impact as before, 
that is g{x) = x and h{x) = x. You may find the optimal trading 
trajectory / for arbitrary A > by first determining /i, using Theorem 



3.2[ next you determine /2 by use of Theorem |3.6[ That is, letting: 

F{f{v)J'{v)) = f'{v)exp{f{v) + f'{v)} 
F\fiv)J'iv)) = FU{v)J-\v)){e'^'^-l) 



Theorem |3.6| states that /2 satisfies the following identity 
-f,{u){2f"{u) + nu) [f'{u) + f"{u)]} 

F\f{v),f{v))dv + 1 



where 
dF^ 



d 



df 



X (^2A 
+2A/'(m 

dF^ 



f2{u) 



dF^ 



d fdF^ 



df dv \ df 



dv = 



du V df 

^f{u)+f'{u) 

X - [{2f"{u) 



f{u){f{u) + nu))] (e" - 1) + (1 + f{u))e-] 
4. Examples 



4.1. Example. The first example we want to study is the case in which 
both the temporary and the permanent impact are linear as in Exam- 



ples |3.3| and |3.7[ The motivation is to compare our results with those 
obtained by Almgren and Chriss (2000). In this example we have set 
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1 and K = 3. The solutions have been numerically calculated 



using Theorems 3^ and 3^ and then plotted in Figures 1 and 2. 

It may be observed that, as one would expect, the solutions are 
not the same and in fact — under the present conditions — our strategy 
dominates Almgren and Chriss. 

4.2. Example. The next example we want to study is the case in 
which the permanent impact is some power less than 1 . Na mely h{x) = 
x^^^, and the permanent is linear. First from Theorem 3.2 we have that 



F{Mu)J[{u)) := sf[{u)exp + {f[{u)f'} 

and /i is the solution to: 

fi - [fi + 3 " (5) " 5^^^^'^' + 5 = °' 

/i(0) = and fi(T) = K. In particular, with T = 1 and i^' = 3 as 
in the previous example we have plotted our result in Figure 3. Note 
that the sublinear impact has increased the speed of execution. Next, 
we computed the Markowitz-optimal open-loop trajectory, with A = 1 
by first computing /2 as described in Theorem |3.6 



Remark 4.1. The previous exercise suggests that if one chooses the 
temporary impact to be sub-linear, the solution — with all the other pa- 
rameters fixed — will always dominate its linear counterpart. On the 
other hand, if the temporary impact is super-linear, the solution will be 
dominated by its linear counterpart. 

A natural question is: What is the correct form of h given our model? 



5. Remarks on Markovian controls 

As pointed out in the introduction, we have only dealt with differen- 
tiable deterministic controls — also known as open loop controls. Fur- 
thermore, our criteria of optimal is in the Markowitz sense. A couple of 
natural and reasonable question arise: how can we study feedback con- 
trols'! How can we optimize with respect to general utility functions? 
Namely, given some utility function U we want to find a trajectory c 
such that 

snpEt,4U{Y)] 

where lA is the set of admissible controls and Y = ^ — KT is the 
execution shortfall. But ^ is in general a very difficult creature to 
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characterize, unless you observe that you may construct a diffusion 
with the following dynamics: 



dXt = ct + 



9{ct) + ^(ct) 



that is equal in distribution to i.e. 

6 = Xt, yt. 



Xt]dt + aXtdBt, = 



Thus 



supEt,,[C/(r)] = supEi,,[C/(XT - KT)] 



and now you may proceed to derive the Hamilton-Jacobi-Bellman equa- 
tion. It is precisely this question, which the authors are investigating 
presently. 
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Figure 1. The graph is plotted with Mathematica. For 
A = the blue line (upper line) represents our optimal 
trading trajectory, the red line is Almgren and Chriss. 
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0.2 0.4 0.6 0.8 1.0 

Figure 2. The graph is plotted with Mathematica. The 
upper red hne is with A = 1, the middle blue line is with 
A = and the lower yellow line is the case of arithmetic 
Brownian motion. 
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0.2 0.4 0.6 0.8 1.0 

Figure 3. The graph is plotted with Mathematica. The 
upper red hne represents the trading trajectory with sub- 
hnear temporary impact. The lower blue line represents 
the trajectory with linear temporary impact. 
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0.2 0.4 0.6 0.8 1.0 

Figure 4. The graph is plotted with Mathematica. The 
upper blue line represents the trading trajectory with 
sub-linear temporary impact and A = 1. The lower red 
line represents the trajectory with sub-linear temporary 
impact and A = 0. 



